Scheme for direct measurement of a general two-qubit Hamiltonian. 



in 
o 
o 

> 

o 

00 
(N 

> 
00 



in ■ 
o ■ 

^ : 

9 L,: 
^— > ■ 



13 



Simon J. Devitt, Jared H. Cole, Lloyd C.L. Hollenberg 1 

1 Centre for Quantum Computer Technology, School of Physics 
University of Melbourne, Victoria 3010, Australia. 
(Dated: February 1, 2008) 

The construction of two-qubit gates appropriate for universal quantum computation is of enormous impor- 
tance to quantum information processing. Building such gates is dependent on accurate knowledge of the 
interaction dynamics between two qubit systems. This letter will present a systematic method for reconstructing 
the full two-qubit interaction Hamiltonian through experimental measures of concurrence. This not only gives a 
convenient method for constructing two qubit quantum gates, but can also be used to experimentally determine 
various Hamiltonian parameters in physical systems. We show explicitly how this method can be employed to 
determine the first and second order spin-orbit corrections to the exchange coupling in quantum dots. 
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The DiVincenzo criterion for quantum computing 1 1 ] em- 
phasises the importance of a universal set of gate operations 
for any physically achievable quantum computer. Two qubit 
operations are an essential component of any universal gate 
library, creating entanglement between qubits. This entangle- 
ment is what gives quantum computers extraordinary power 
over classical computational devices. The design of appro- 
priate two qubit operations requires accurate knowledge of 
the interaction dynamics between qubits, specifically the two 
qubit interaction Hamiltonian. Accurate knowledge of this 
Hamiltonian is especially important for solid state devices, 
where qubit fabrication could easily lead to interaction dy- 
namics which vary from qubit to qubit. 

While in theory, many solid state systems rely on a purely 
Heisenberg type Hamiltonian [^blcBi anisotropic terms can 
be present when taking into account higher order effects. For 
example, spin-orbit coupling in quantum dots j^, 0, 0] where 
the Hamiltonian takes the time-/ndependant form, 



H = J[Si ■S 2 +P-(S 1 x S 2 ) + Sirs 2 ] 



(1) 



First order corrections are represented by j3, known as the 



ited by 

Dzyaloshinski-Moriya (DM) vector |8, 9] and T is a rank-3 
symmetric tensor representing second order corrections. Al- 
though theoretical estimates on these corrective terms have 
been made |7], and measurements of the DM corrections 
made in spin glass systems 1 10], there has been no experimen- 
tal measurement of either (3 or T for isolated coupled dots. 

In the context of quantum information processing, The sys- 
tem Hamiltonian is used to design quantum gate operations. 
A large amount of work has already been completed regard- 
ing the accurate construction of two qubit gates 

113 113 

Il4lll5ll . this work can be split into two broad categories. The 
first assumes that the quantum computing system exhibits a 
discrete well defined two qubit interaction. This occurs for ar- 
chitectures such as linear optics systems Ill6 lll7ll . with CNOT 
gates experimentally demonstrated 1181 [1911 . These discrete 
two qubit operations are traditionally characterised using state 
and process tomography II19l |20j |2ll 12211 where the interac- 
tion Hamiltonian is initially assumed on well founded theo- 
retical and/or experimental grounds. Arbitrary gates can then 



be constructed directly from this fundamental discrete inter- 
action JT4IH5I1 . While tomography has its advantages and is 
currently the method used by experimentalists, it cannot di- 
rectly extract the Hamiltonian and becomes difficult when the 
gate is not ideal and a good theoretical model is unavailable. 

The second category is exemplified by solid state architec- 
tures. Here two qubit interactions exhibit time dependent con- 
trol where unitary gates, U(t), can be specified by the interac- 
tion Hamiltonian H, U(t) = exp(iHt). State and process to- 
mography can be performed, but in general a discrete gate has 
to be designed beforehand that can produce entangled states 
(for example a CNOT). However, it is much more efficient to 
design interactions directly, using the Hamiltonian [11]. 

Methods for both single and two qubit Hamiltonian char- 
acterisation has recently been developed I23ll24l|25)]. These 
methods involve mapping the system evolution 1 23j l24ll (sin- 
gle qubit characterisation) or the entanglement of the system 
112 511 (two qubit characterisation) over time. However, pre- 
vious work on two qubit characterisation restricted the inter- 
action Hamiltonian to a Heisenberg form, H = c\XX + 
c 2 YY + c 3 ZZ, where XX = a x ® a x , YY = a y ® a y 
and ZZ = u z ® a z . In this letter, we present a more general 
method for characterisation. We will show how mapping the 
entanglement of the two qubit system gives enough informa- 
tion to not only determine the entangling properties of a given 
interaction, but to fully reconstruct the Hamiltonian. We ap- 
ply this method to characterise a trial Hamiltonian of the form 
shown in Eq. ^and therefore determine both (3 and T. 

The most general form of a fully non local (FNL) two qubit 
Hamiltonian is H = Y^i j=i dijO~i®ojj where {ai, 02, 03} = 
{a x , Uy, a z } and dij G M. Ref. [11] shows that for any 
dij there exists single qubit operators K G SU (2) ® SU (2) 
such that U(t) = exp(iHt) = K^exp(iH a t)K, where H a = 
Cl XX + c 2 YY + c 3 ZZ and [c u c 2 ,c 3 ] G E. If H a and if can 
be experimentally determined, then the original Hamiltonian 
can be reconstructed, H = K^H a K. 

Determining H a and K can be done by examining the en- 
tanglement produced by the operator U(t) on a known sepa- 
rable input state \ip). The entanglement measure of the state 
\4>) = U\tjj), and hence U can be quantified by the func- 
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tion C 2 — \((f)*\YY\(j))\ 2 , denoted the squared concurrence 
lEill . C 2 varies continuously between 0, for a product state 
and 1, for a maximally entangled Bell state. Experimental 
techniques have been developed to measure this function for 
a given two qubit state that requires measurements in the ZZ 
and XZ basis lE5l Eill . Our method for using C 2 to recon- 
struct the Hamiltonian will assume that single qubit gates have 
been characterised to arbitrary accuracy using the method de- 
tailed in 12411 . hence we have full single qubit rotation control 
and that measurements can be performed in these two bases. 

Consider an unknown Hamiltonian H = Xh=i <kj Pi ® o~j- 
U(t) — exp(iHt) can be expressed according to 111 ill as 
U(t) — K^exp(iH a t)K . Given a known separable input 
state C 2 for the resultant state U(t)\4>o) can be writ- 

ten as C 2 = \(Q^K\ipo}) T F 2 (Q^K\ip }\ 2 , where F = 
di ag {e iAlt ,e a2t ,e a3 *,e a4t } and {A x , A 2 , A 3 , A 4 } = {ci - 
c 2 + c 3 , ci + c 2 - c 3 , -ci - C2 - c 3 , -ci + c 2 + c 3 }. The 
matrix Q transforms computational states to Bell states full . 

The added complexity in characterising FNL Hamiltonians 
compared with the non-isotropic Heisenberg Hamiltonian de- 
tailed in l25ll manifests itself in the K matrix present in the 
expression for C 2 . In the Heisenberg case, K = II, where II 
is the 4x4 identity matrix. For FNL Hamiltonians, we have 
no knowledge of H, hence we have no knowledge of K. The 
initial known product state |^o) is therefore rotated by K be- 
fore it is acted upon by the interaction term exp(iH a t). Unlike 
characterising Heisenberg Hamiltonians, \ipo) cannot be cho- 
sen to isolate [ci, c 2 , c 3 ] directly. 

Consider a general 4x4 matrix K e SU(2) ® SU(2). 
It is known that the group SU(2) (g> SU(2) is isomorphic to 
50(4) |11], in fact SO (A) = Q^SU{2) <g> SU{2)Q. Hence 
Kq = KQ £ 50(4) and all elements in Kq are real. Let 

q%) = K Q Q^\ib Q ) = h\9+)+ia\^-)+h\9+)+u\9-), 

be the rotated input product state in the Bell basis. In this 

case, the concurrence can be rewritten as C 2 — \e 2tXlt l 2 + 
e aiA,t i a + e aiA 8 t f a + e aiA4t i a|a i 

Experimentally, two qubits are initialised in a known prod- 
uct state \ipo) an d C 2 mapped out as a function of t. Once C 2 
has been mapped the data will have the following form, 

^=x>t+e \m 2 M {uiit+iii) > (2) 

i— 1 i>t3~>i 

where u>ij = 2(A, — Aj) and 7^ = 2(arg(Zj) — arg(Zj-)). The 
Fourier spectrum of the time series, C 2 , is taken and the power 
density spectrum |g(w)| 2 is plotted as a function of uj. 

We will restrict this analysis to the case where H decom- 
poses to a form where c\ 7^ c 2 7^ c 3 7^ 0. This represents the 
most straightforward case for characterisation and arguably 
the most likely given a random Hamiltonian. Other cases, 
i.e. when c\ — c 2 = c 3 , c 3 = etc, can still be charac- 
terised, however the method is slightly more involved. The 
power density spectrum, \g(u)\ 2 , has the form, 

(2tt)| 5 H| 2 = (j2 \lA 5(u) + £ \h\%fS^, (3) 

i— 1 hj^i 



with <5 Wi . = [<5(w + u>ij) + 6(u> — Ldij)]. All the information 
on the rotated input state KqQ'\iPq) is contained within the 
height of each peak, while all the information relating to the 
factors [c% , c 2 , c 3 ] are related to the location of each peak in 
the frequency spectrum. 

Extracting the actual matrix Kq and [ci,c 2 ,c 3 ] from the 
power density spectrum of Eq. [3] is reasonably straightfor- 
ward when c\ 7^ c 2 7^ c 3 7^ 0. Identifying each of the six 
peaks in the frequency spectrum relies on the principal of lo- 
cal equivalence II ill . Decomposing any H into H a defines 
a local equivalence class [ci, c 2 , c 3 ] since any two operators 
t/i and f/ 2 which are locally equivalent differ through single 
qubit rotations. The values of [ci, c 2 , c 3 ] are therefore locally 
periodic. For example the single qubit operators ±iX (g) X 
permutes [ci,c 2 ,c 3 ] — > [ci ±7r/2, c 2 , c 3 ], and similarly for the 
local operators ±iF ® Y and ±iZ ® Z. This locally equiv- 
alent periodicity allows the restriction of [ci, c 2 , c 3 ] such that 
ci > c 2 > c 3 > 0. Enforcing this condition allows the identi- 
fication of all 6 separate peaks within the Fourier spectrum. 

From the definition of Uij we find, cji 2 = 4(c 2 — c 3 ), 
W13 = 4(ci + c 3 ), W14 = 4(ci - c 2 ), cj 23 = 4(ci + c 2 ), 
^24 = 4(ci — c 3 ) andw 3 4 = 4(c 2 + c 3 ) (overall negative signs 
are omitted since Eq. [5]is symmetric in coy — > —u>ij). Using 
the restriction of ci > c 2 > c 3 > and ci 7^ c 2 7^ c 3 7^ 0, 
the two largest frequencies for all c, are w 23 > u>i 3 , which are 
separated by u>i2 = 4(c 2 — c 3 ). w 2 4 > W14 is also separated 
by lo\2- All four frequency peaks, w 23 , u>\%, uj\ 4 and w 2 4 are 
cantered on the value Ac\. The final peak is cj 3 4, which with 
uj\2 is symmetrically spaced about 4c 2 < 4ci, (see Fig. [Q. 

The coefficients [ci, c 2 , c 3 ], which can be used to define H a 
can be determined, however to reconstruct H, K such that 
H = K'H a K also needs to be found. This can be done 
through the peak heights in the Fourier spectrum. Consider 
the matrix Kq = Q^KQ S 50(4). All elements in Kq are 
real and six independent parameters are required to define an 
arbitrary 50(4) matrix. In the case where ci 7^ c 2 7^ c 3 7^ 
parameterisation of Kq in order to determine these six pa- 
rameters isn't required. Instead, all 16 elements of Kq can be 
found directly. Consider the following three separable states 
(neglecting normalisation), |i/>o)i = (i|0) + |l))®(|0) + i|l)), 
l^o)2 = (|0) + |1» ® |1), |^o)3 = |1) ® (|0) + If 
we specify the matrix Kq through the 16 matrix elements 
{ai, ..,,a\§}, then the squared modulus of the rotated Bell 
state coefficients, \lf(k)\, for each of the three input states, k, 
are, 2|/ 2 (1)| = (a 2 + a 2 ), 4|/ 2 (2)| = (ai + a 3 ) 2 + (a 2 + a 4 ) 2 
and4|Z 2 (3)| = (ai +a 3 ) 2 + (a 2 — a4) 2 . Similar expressions 
are obtained for |/|(fe)|, |^ 3 (^) and |^ 2 (fc)| using the appropri- 
ate rows from the matrix, Kq. 

Examining Eq. [3] for the power density spectrum of C 2 , 
the height of each of the respective peaks are, Pij(k) = 
Zi(fc)| 4 |(j(fc)| 4 . Since \h(k)\ > 0, we can determine the ab- 
solute value of the rotated input state (in the Bell basis) as, 

I/, = / P12P13 \ = ( P12P14 \ = f PuPiz \ 

1 \ P23 / V P24 / V P34 / 
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/ \ 1/8 / \ 1/8 / \ 1/8 

/ P12P23 \ _ / P12P24 \ _ / P24P23 \ 

V P13 J \ P14 / V P34 / 

/ \ 1/8 / \ 1/8 / \ 1/8 

/ P13P23 \ / P13P34 \ / P23P34 \ 



V Pl2 J V Pl4 J \ P24 / 

/ \ 1/8 / \ 1/8 / \ 1/8 

1^1 ( P1AP24 \ = /P24P34\ = / P14P34 \ 

4 \ P12 / \ P23 ) V Pl3 J 

(4) 

dropping the index k. Invoking the group conditions of 
50(4), specifically, that the 4 rows of Kq forms an orthonor- 
mal set, leads to the following, 



/ , 1 (1) = |/ 1 (1)| 2 -|Z 1 (2)| 2 + |Z 1 (3)| 2 = 
^(l)-|^i(l)| 2 + |/i(2)| 2 -|Z 1 (3)| 2 = 

Mi) = -IMi)l 2 + IM2)l 2 + IM3)l 2 

M4 (l) = l-|Zi(l)| 2 -|Zi(2)| 2 - 



~(a 3 -ai) 2 , 
~(a 3 +ai) 2 , 

= ^(&2 - A4) 2 , 



(3)| 2 = -(o 2 + a 4 ) 2 . 



(5) 

Similar equations can be used to determine the other three 
rows of Kq. Since n represents the square of the matrix 
elements, there is a certain amount of ambiguity for each row 
of Kq. Eq. |5] give rise to solutions in pairs, (01,03) and 
(02,04) which implies that various matrices can be formed 
for Kq which vary through permutations and/or sign flips of 
these pairs. This ambiguity leads to 2 16 matrices for Kq. 

Noting that Kq € SO (4) eliminates a large number of 
possibilities. Each possible matrix satisfying Eq. [5] can be 
numerically constructed and checked that KqKq = II and 
Det(ifg) = 1. However, even in the case where all the 
elements of Kq are distinct, there still exists a number of 
matrices satisfying the group conditions of 5*0(4). Each of 
these matrices, Kq, combined with the already determined 
H a , will lead to a set of locally equivalent Hamiltonians, 
{Hi}, i.e. Hamiltonians with the same values for [ci, c 2 , C3]. 

In order to narrow down the possible Hamiltonians, we 
analytically construct all the operators, {Ui(t) = exp(ii^i)} 
for each locally equivalent Hamiltonian. For each Ui(t), C 2 
can be calculated for all three input states \ipo)k at some 
arbitrary value of t. Comparing the experimental data to the 
analytic calculations of O 2 will reduce the possible lists of 
Hamiltonians to a smaller subset that have precisely 

the same analytical form for C 2 , for all input states |i/ , o)fc- 

To isolate the correct Hamiltonian from this finite set, 
{Hi/}, analytically construct the function C 2 for some arbi- 
trary input state j^o) that produces distinct functional forms 
(for C 2 ) for each XJy (t). Once the appropriate input state has 
been analytically determined, experimentally measure C 2 
for some single arbitrary value of t. The experimental and 
analytical value of C 2 can be compared for each £7^ and the 
correct Hi> determined. Once this procedure has been com- 
pleted, only the matrices, Kq, that produce the same target 
Hamiltonian will remain, satisfying all the experimental data. 



Although the characterisation process is essential to the ef- 
ficient construction of arbitrary two qubit gates, it can also be 
used to experimentally determine other theoretical parameters 
governing a system. The following is a specific example to 
highlight the methodology. Consider the exchange interaction 
between spin systems in quantum dots. As mentioned in 
the introduction, spin-obit coupling introduces anisotropic 
corrections to the otherwise ideal isotropic exchange Hamil- 
tonian, H = JSi.S2- When these correction are included, 
the time independent Hamiltonian describing the two qubit 
coupling takes the form in Eq. \l\ Although gate design 
schemes have been developed to overcome these anisotropic 
terms |6], characterisation of this coupling is possible and 
allows us to determine both (3 and T up to a factor of J. 

We will consider the trial Hamiltonian where 
J = 1, = {0.01,0.005,0.02} and T = 
{0.003, 0.005, 0.001, 0.0015, 0.0024, 0.0009} (represent- 
ing the 6 independent parameters specifying a 3 x 3, real, 
symmetric matrix). These values are motivated by the fact 
that J = 1 leads to extracting (3 and V exactly, f3 is chosen 
randomly but is consistent with the analysis of Kavokin Q] 
for quantum dots in GaAs. T is roughly an order of magnitude 
smaller since it is a second order correction to H. 

Fig. [0 shows the power density spectrum obtained using 
the input state ^0)2 = (|0) + |1))|1). From the spec- 
trum we can identify the relevant peaks as W23 = 8.039, 
W13 = 8.025, W12 = 0.012, cj 24 = 0.028, w 14 = 0.016 
and CJ34 = 8.011. These peak frequencies lead to 
[ci,c 2 ,c 3 ] = [1.006,1.003,0.999]. The respective peak 
heights and values for (pi,fi 2 , ^3,^4) are given in tableU 

After constructing all possible combinations of the above 
values for {ax, ai6}> 64 matrices are found to satisfy the 
group conditions of 50(4). These matrices and H a are used 
to find 16 distinct Hamiltonians H = K^ H a K. Eq. |6]shows 
4 Hamiltonians in terms of the spin-orbit corrections (3 and V. 



H 1 : /3= {-0.0015, 0.005, -0.009} 

T = {-2, 0.005, -2, -0.01, 0.002, -0.002} 

H 2 : 13 = {0.01,0.005,0.02} 

T = {0.003,0.005,0.001,0.0015,0.0024,0.0009} 

H 3 : j3= {0.009, 0.005, 0.0015} 

T = {-2, 0.005, -2, 0.02, 0.002, 0.01} 

Hi : j3= {0.02,0.005,0.01} 

T = {0.001, 0.005, 0.003, 0.0009, 0.002, 0.0015} 

(6) 

When parameterising Eq. [6]with respect to (3 and T, we have 
explicitly used J = 1, this is NOT a requirement to charac- 
terise the Hamiltonian using this method. However, since the 
general form of the Hamiltonian has J as a constant multi- 
plicative factor, extracting the exact values for (3 and T re- 
quires J to be known or assumed. 

From the 16 possible Hamiltonians generated, clearly the 
Kq matrices generating H 2 are the target matrices, how- 
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FIG. 1: Power Density spectrum for the |i/»o)2 input state. Note that 
the respective symmetries discussed earlier. The values for each peak 
height are shown in table Q. 
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P23 


P24 
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6.81 


2.32 








13.35 




I0O>2 


.0467 


5.68 


11.01 


.0197 


.0383 


4.66 




I0O>3 


.0473 


5.91 


10.84 


.0201 


.0367 


4.59 




Ml(l) 


Mi) 


Ms(l) 


M4(l) 


Mi (2) 


M2(2) 


M3(2) 


M4(2) 


21.21 


20.46 


56.01 


2.33 


.0041 


.021 


4.49 


95.48 


Mi (3) 


M2(3) 


Ms(3) 


/i 4 (3) 


Mi (4) 


Ma (4) 


M3(4) 


M4(4) 


50.24 


49.74 


.0110 


.0012 


28.55 


29.78 


39.49 


2.19 



TABLE I: Relevant data from the three power density spectra 
in Fig. Q Values for all the respective peak heights and 
{/Lti (i) , 1x2 {%) , M3 {i) , M4 (i) } can be used to find i^Q . 



ever we need to filter out the other 15 possible Hamiltoni- 
ans. To do this, analytically construct the unitary operators 
Ui = exp(iHi), for i £ [1, 16] and calculate C 2 for all 
i £ [1, 16] and k £ [1,3]. The calculated values are then com- 
pared with the values of C 2 measured experimentally at t = 1 
for all three input states. Table UTI shows the results for the 4 
Hamiltonians in Eq. [6] Simulated values of C 2 immediately 
eliminate Hi and H% as possible Hamiltonians. Repeating for 
the other 12 Hamiltonians shows that the analytical values for 
C 2 match only for H2 and H4. Here, the analytical forms for 
C 2 are identical for the three input states \4>o)k- 

To isolating H2, we analytically construct C 2 and find that 
the input state | tp ) 4 = 1 1) ( 1 0} — i\ 1) ) generates different func- 
tions for H2 and H4. One more experimental measurement of 
C 2 is taken for this input state at t — 0.5, giving C 2 = 0.200. 
Analytically, C 2 (H 2 ) = 0.200 and C 2 (iJ 4 ) = 0.196. This fi- 



Input 


Ent(£/i) 


Ent(t7 2 ) 


EBt(C/ 3 ) 


Ent(£/ 4 ) 


Ent([/(1)) 


l</>o)i 


.581 


.581 


.581 


.581 


.581 


I0o)2 


.147 


.147 


.142 


.147 


.147 


I0O>3 


.144 


.150 


.150 


.150 


.150 



TABLE II: Calculations of C 2 for all three input states using the 4 
unitary operators experimentally obtained fort = 1. The last column 
represents the data for the unknown unitary operator at t = 1, U(l). 



nal measurement has discriminated between H2 and H4, leav- 
ing H2 to satisfy all the experimental measures of C 2 . 

In conclusion, we have presented a systematic method to 
experimentally determine the interaction Hamiltonian of a 
two-qubit system that requires initialisation in a minimum of 
three product states and measurement in two separate bases. 
We have demonstrated how this method not only allows for 
accurate characterisation of system Hamiltonians required for 
constructing two-qubit quantum gates, but also how charac- 
terisation can lead to the determination of other parameters 
of interacting systems. Specifically, we showed how the first 
and second order anisotropic corrections, in quantum dots can 
be experimentally determined. A comprehensive error analy- 
sis of this scheme under a fixed number of measurements and 
how mis-characterisation propagates through to effective sys- 
tematic gate errors is still required. 
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